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Abstract. Wc present a detailed study of the spectral and temporal properties 
of the X-ray and optical emission of GRB050713a up to 0.5 day after the main 
GRB event. The X-ray light curve exhibits large amplitude variations with several 
rebrightenings superposed on the underlying three-segment broken powerlaw that 
is often seen in Swift GRBs. Our time-resolved spectral analysis supports the 
interpretation of a long-lived central engine, with rebrightenings consistent with 
energy injection in refreshed shocks as slower shells generated in the central engine 
prompt phase catch up with the afterglow shock at later times. Our sparsely- 
sampled light curve of the optical afterglow can be fitted with a single power law 
without large flares. The optical decay index appears flatter than the X-ray one, 
especially at later times. 

Key words, gamma ray:bursts- gamma ray: individual GRB 050713A 
1. Introduction 

A commonly accepted wisdom in the pre-Swift era was that the optical and X-ray after- 
glows of Gamma Ray Bursts (GRBs) showed a smooth power law decay with time after 
the burst (e.g. Laursen & Stanek 2003). This behavior was found to be consistent with 
the fireball model, where the afterglow flux is produced when a relativistic blast wave 
propagating into an external medium. Under the assumption of a spherical fireball and a 
uniform medium. Sari, Piran & Narayan (1998) showed that the dependence of the flux 
on frequency and time can be represented by several power law segments, oc u~^t~°'. 
Before Swift only a few GRBs showing deviations from the smooth power law light curve 
were known, see e.g. the case of GRB 021004 (Bersier et al. 2002; Matheson et al. 2002). 
Its lightcurve was densely sampled in the optical, allowing detailed modeling. The bumps 
were interpreted as due to overdensities in the interstellar medium in which the after- 
glow is produced (Lazzati et al. 2002; Nakar et al. 2003; Heyl & Perna 2003). On the 
other hand, Bjornsson et al. (2004), Nakar et al. (2003) and de Ugarte Postigo et al. 
(2005) modeled these fluctuations as due to several energy injection episodes. Several 
re-brightenings were observed also in the optical light curve of GRB 030329 and were 
interpreted as due to refreshed shocks (Granot, Nakar & Piran 2003, Huang et al. 2006). 

This simple picture is now changing since the advent of Swift. Bright X-ray flares 
have been recently observed by Swift in almost half of its detected GRBs (Gehrels et al. 
2005; Burrows et al. 2005b; Falcone et al. 2006; Nousek et al. 2005; O'Brien et al. 2006). 
While some bursts show one distinct flare, like GRB 050406 (Romano et al. 2006), other 
events like GRB 050502B and GRB 050713A show several flares (Burrows et al. 2005; 
O' Brien et al. 2006; Falcone et al. 2006). One of the main current goals of the GRB 
community is to understand the origin of this newly observed light curve behaviour. Since 
flares are likely to trace the activity of the internal engine, they can help us gain a more 
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comprehensive view of the physical processes governing the early phases of GRB activity. 
Flares have been observed in the light curves of both long and short GRBs. For the long 
bursts, King et al. (2005) proposed a model in which the flares could be produced from the 
fragmentation of the collapsing stellar core in a modified hypernova scenario. For the case 
of short GRBs, MacFadyen et al. (2005) suggested that the flares could be the result of the 
interaction between the GRB outflow and a non-stellar companion. More recently, Perna 
et al. (2006) have analyzed the observational properties of flares in both long and short 
bursts, and suggested a common scenario in which flares are powered by the late-time 
accretion of fragments of material produced in the gravitationally unstable outer parts 
of the hyperaccreting accretion disk. Other mechanisms that could produce flares are of 
magnetic origin (Gao et al. 2005; Proga & Zhang 2006). In this paper we concentrate on 
one particular burst that presents flaring activity, GRB 050713A, and perform a detailed 
spectral and timing analysis of its main flare with the goal of constraining the physical 
mechanisms that can be responsible for its production. 

The Swift BAT localized this burst on 2005, July 13.1866 UT to a 3' radius error 
circle (Falcone et al. 2005). The BAT light curve is characterized by a main event lasting 
~ 13sec that drops by a factor of 100, followed by two rebrightenings starting ~ 53 s and 
~ 110 s after a time to s, that corresponds to the time at which BAT started to detect 
the burst. The spectrum of the main event could be fitted with a power law with energy 
index 0.58 ± 0.07, yielding a fluence of (9.1 ± 0.6) x 10"*^ erg cm'^ (15-350 keV). The 
peak flux in 1 s window time is 6.0 ± 0.4 ph cm~^ s~^ (Golenetskii et al. 2005). Swift 
slewed promptly toward the position of the GRB and XRT started to observe this event 
just 70 s after the trigger. The XRT light curve was dominated by a major rebrigthening 
event 117 s after to, nearly coincident with the second rebrightening detected by BAT. 
Both the 15-350 keV and the 0.5-10 keV fluxes varied up and down by a factor of ~ 10 
in about 40 s. A second XRT rebrightening event occurred about 186 s after the to, and 
other smaller amplitude flares were seen at later times (around ~ 10* s). In this paper 
we focus on the first major XRT rebrightening event through both a temporal and a 
time-resolved spectral analysis. The statistics of the spectra of the second flare are not 
good enough to allow a similar detailed analysis. Furthermore, we use XMM-Newton 
data to better constrain the late time afterlow decay. 

We show that the rebrightening in the early-time afterglow light curve is consistent 
with energy "injections" , probably due to later shells that catch up with the main shock 
at later times. 

The redshift of GRB 050713A is not known. The host galaxy is not detected, and the 
limit on its magnitude is quite shallow (i?^ 23), due to the high background caused by a 
nearby (~ 1 arcmin), bright star. 
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2. Observations and data reduction 



2.1. Swift XRT observations 

XRT started observing GRB050713A at 2005-07-13, 04:30:14.9, just 70 s after the BAT 
trigger. The observation lasted mitil 2005-07-13, 08:50:27 (UT), for a total of - 1700 s net 
integration time. Within this observation, XRT observed the GRB in three consecutive 
orbits. 

The data were reduced using the Swift Software (v. 2.0) and in particular the XRT 
software developed at the ASDC and HEASARC (Capalbi et al. 2005^). Standard screen- 
ing criteria were adopted to reject "bad" events following Capalbi et al. (2005). In partic- 
ular, the constraint on the angular distance between the source position and the satellite 
position (which must be less than 0.08 deg, Capalbi et al. 2005), led us to exclude the 
first part of the first orbit and 80% of the last orbit. The GRB was observed by the XRT 
in two observing modes. During the first satellite orbit the GRB was mostly observed in 
windowed timing mode (WT mode, providing ID imaging). We report here the WT data 
only for this orbit. During the second and third orbits the GRB was mostly observed 
in photon coimting mode (PC mode, providing the usual 2D imaging). We report here 
the PC data only for these orbits. We used a — 12 grade selection for the PC mode 
and — 2 for the WT mode. We excluded all the events with energy below 0.3 keV, to 
minimize the background due to the bright Earth limb and and because the calibration 
of the data below 0.3 keV are still uncertain. The intensity of the source was high enough 
to cause significant pileup in the PC mode. In order to avoid this pileup we extracted 
counts from an amnilus with inner radius of 6 pixels and obiter radius of 20 pixels (14 and 
47 arcsec, respectively). We then corrected the observed count rate for the fraction of the 
XRT Point Spread Function (PSF) lying outside the extraction region. The correction 
was equal to a factor of 3.33 (i.e. 30% of the PSF lying inside the annulus). Data in 
WT mode were not affected by pileup and therefore we extracted counts from a circular 
region of 20 pixels radius (47 arcsec). Physical ancillary response files were generated 
with the task xrtmkarf to account for the different extraction regions. For the spectral 
fits we used the latest redistribution matrices (version 7). 

2.2. XMM-Newton observations 

XMM-Newton (Jansen et al. 2001) observed this GRB for about 30 ks starting from 2005- 
07-13 10:18 UT, about six hours after the trigger, observation ID 0164571001 (Loiseau 
et al. 2005). The data were processed using the XMM-Newton Science Analysis Survey 
(SAS) v. 6. 1.0^. We used the raw event files (i.e. the observation data files, ODF), which 

^ http:/ /heasarc. gsfc.nasa.gov/docs/swift/analysis/xrt_swguide_vl_2.pdf 
^ http://xmm.vilspa.esa.es/external/xmm_ sw_cal/sas_frame.shtml 
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were linearized with the XMM-SAS pipehnes, epchain and emchain for the PN (Striider 
et aL 2001) and MOS (Turner et al. 2001) cameras respectively. Events spread at most 
in two contiguous pixels for the PN (i.e. grade=0 4) and in four contiguous pixels for the 
MOS (i.e, grade=0 12) have been selected. Event files were cleaned from bad pixels (hot 
pixels, events out of the field of view, etc.). In order to remove periods of high background, 
we analyzed the light curves of the counts from the entire EPIC PN and MOS CCDs 
at energies higher than 10 keV, where the X-ray sources contribution is negligible. We 
rejected time intervals, in which this comit rate was higher than 10 counts s^^ and 1.5 
counts s~^ for the PN and MOS cameras respectively. This corresponds to rejecting 
27% (15%) of the PN (MOS) on source time. The source counts were extracted from 
a circular region of 47 arcsec radius (12 pixels). The background coimts were extracted 
from the nearest source free region. The response and ancillary files were generated by 
the XMM-SAS tasks, rmfgen and arfgen respectively. 

3. Optical observations 

UVOT started observing the field at the same time of XRT (75 seconds after the BAT 
trigger) but did not detect any new source in the XRT error circle, with a 3a upper 
limit of 18. Shortly after the receipt of the alert, we started observing the GRB 
location in order to look for an optical counterpart. Observations were conducted using 
the Liverpool Telescope (LT), the Telescopio Nazionale Galileo (TNG), and the Nordic 
Optical Telescope (NOT), all located in the Canary Islands. A bright object was discov- 
ered inside the XRT error circle by inspecting the TNG 7-band images, taken 50 min 
after the GRB, and was proposed as the GRB afterglow (Malesani et al. 2005). Its coor- 
dinates were aj2ooo = 21 : 22 : 09.57, (5j2ooo = +77 : 04 : 29.4 (0.15" uncertainty). Hearty 
et al. (2005) subsequently reported variability for this source, confirming that it was the 
GRB counterpart. The robotic 2-m Liverpool telescope imaged the GRB field starting 
from 150 seconds after the burst, detecting the afterglow in the r' band (Monfardini et 
al. 2005). Further imaging was secured at the NOT in the R band (Jul 13; Guzly et al. 
2005) and at the TNG in the / band (Jul 14). Only an upper limit could be obtained at 
the Observatorio de la Sierra Nevada (OSN) during the night of Jul 13. Unfortunately, 
the GRB position was close (ss 1') to a bright foreground star {R ~ 6.5), which caused 
a bright sky level, so that no further detection was possible in the following nights. 

The TNG and NOT data analysis and reduction were carried out following the stan- 
dard procedures. The LT data analysis was performed using the dedicated LT pipeline 
(Guidorzi et al. 2006). Magnitudes were computed using the SExtractor (Bertin & 
Arnouts 1996) and Gaia software packages. Photometric calibration was obtained by 
observing standard fields at the TNG (7 band) and OSN {R band). The zeropoint, how- 
ever, was computed at only one epoch, so it may suffer from systematic uncertainty. LT 
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Fig. 1. The BAT XRT and XMM-Newton PN light curve of GRB 050713 A and its 
afterglow (solid line and points with errors). The BAT 15-300 keV count rate, the XRT 
and PN 0.5-10 keV count rates are converted to a flux at 1 keV assuming a power law 
spectrum with average spectral indices /3 = 0.5 (BAT) and /? = 1 (XRT, XMM). Filled 
circles are the optical observations performed with the RAPTOR-S, Liverpool 2-m, TNG 
and NOT telescopes. The dashed line is the best fit power law decay of the optical data. 



data were calibrated via Landolt standards with SDSS calibration (Smith et al. 2002). 
We combined all our data, and included the early RAPTOR i?-band detection (Wren et 
al. 2005). A log of all optical observations is given in Table 1. 

Figuren]shows the light curve of the optical flux at TOOOA obtained by converting the 
observed R, r' and I magnitudes into flux at TOOOAassuming a power law with spectral 
index of -1. Considering a power law flux decay F{t) oc we found a = 0.67 ± 0.05. 
We note that, to within the limitations of the sampling, the optical light curve appears 
fairly smooth, with no signs of the large fluctuations seen in the X-ray region. However, 
the coverage is quite scarse, so that no strong conclusion can be drawn. 
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Table 1. Log of Optical observations 
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4. Results and discussion 

4.1. Light curves 

Figure n shows the BAT, XRT and XMM-Newton Ught curve of GRB 050713 A and its 
afterglow. The power law slope of the XRT decay starting about 200 s after the trigger 
was 0.67, very similar to the optical decay starting 150 s after the trigger (filled circles 
in figure^. The power law decay during the XMM observation is 1.5. Thus, the light 
curve is broadly consistent with that of many Swift afterglows (Chincarini et al. 2005, 
Tagliaferri et al. 2005, Nousek et al. 2005), with a steep decay followed by a shallower 
decay and finally a transition to a standard a — 1.5 decay. The break between the 
shallow and steep decay should occur beween ~ 5000 s and ~ 20.000 s otherwise the 
extrapolation of the flux, based on the XMM decay index, would largely exceed the flux 
actually detected by XRT in PC mode. 

The BAT data are characterized by two rebrightenings ^ 53 s and ^ 110 s after tg. 
The early XRT-WT data are characterized by prominent flares ~ 128 and 188 seconds 
after tg. We note that the first XRT flare overlaps with the second BAT flare, and that 
the peak in the hard X-ray band is shifted of 10 s respect to the soft one. An additional 
flare is clearly visible in the XRT PC data. In the following, we discuss in more detail 
these features. Figure [3 shows a zoom of the light curve of the flrst flare in three energy 
bands: 0.8-1.4, 2.2-4 and 4-10 keV. Three features are apparent from this flgure: (a) 
the rise time looks similar in all energy bands; (b) the decay of the flrst flare is sharper 
at higher energies; (c) the peak of the softer X-ray light curve is shifted by about 10 s 
with respect to that of the harder X-ray light curve. To make these statements more 
quantitative, we computed the power law rise indices Ur from 110 to 118 s after to and 
the power law decay indices ad from 124 to 156 s after to, for the light curves in the 



8 



Guetta et al.: GRB 050713A: the case of refreshed shocks 




no 120 130 140 150 160 

Time (s) 



Fig. 2. Light curves of the first flare in three energy bands. 



0.3-0.8, 0.8-1.4, 1.4-2.2, 2.2-4 and 4-10 keV bands. We used bins of 4 s for the first 
band and bins of 2 s for the other 4 bands. The rise time indices were all consistent 
with the value 0.6. Conversely, figure 01 shows the best fit decay indices as a function of 
the energy. We note that the power law rise and decay indices depend strongly on the 
initial counting time, Tq. Our choice was to set Tq at the beginning of the flare rise for 
ar and at the beginning of the decay for ad, respectively. Setting Tq at the GRB trigger 
time would produce much steeper indices (e.g. ad > 3), possibly requiring a different 
physical interpretation of the event, such as late internal shocks (Zhang et al. 2005). The 
underlying decay of the X-ray light curve between 70 and 200 seconds after the main 
GRB events is not well defined in this case, but it appears roughly consistent with the 
a ^ 3 decay of Zhang et al. (2005). 

We also computed cross-correlation functions between the light curves in the four 
hardest bands and in the 10-30 keV band with respect to the 0.3-0.8 keV band. We used 
the first 128 s of WT data for this analysis, which includes the two prominent peaks, in 
4 s bins. We fitted the peak of the resulting cross-correlation functions with a Gaussian 
to estimate time lags. Figure 0] shows the time lag as a function of the energy. We have 
also tried to compute the cross-correlation functions after subtracting a power-law trend, 
computed using the first 20 s of WT data and extrapolating the best fit to the rest of 
the WT data. The time lags obtained subtracting or not subtracting the power law trend 
are fully consistent one with the other. 



Guetta et al.: GRB 050713A: the case of refreshed shocks 



- 



-0.5 



-1 



-1.5 



T — I — I — I — r 



T — I — I — r 



_i — — — — — — — 



p = 2 



p = 2.5 



-£3 ' 



J I I I L 



J I I L 



10 



Energy (keV 



Fig. 3. Power law decay indices of the light curves in the 0.3-0.8, 0.8-1.4, 1.4-2.2, 2.2- 
4, and 4-10 keV bands, starting from 108 s after to, as a function of the energy. The 
prediction of the synchrotron model is represented by the dashed lines considering a 
broken power law for two different values of p. 



4.2. Time-resolved spectroscopy 

The analysis of the WT light curves of GRB 050713A revealed complex spectral varia- 
tions. To investigate the nature of these variations we performed a time resolved spectral 



10 



Guetta et al.: GRB 050713A: the case of refreshed shocks 




1 10 

Energy (kcV) 



Fig. 4. The time-lag computed from the hght curves of the first 128 s of XRT-WT 
observations, as a function of the energy. The two sohd curves show the expected time 
lags for Cs"^ £^52^71]^^ = 10^^ (upper curve) and 2 x 10^^ (lower curve) of the standard 



Sari & Piran (1999) afterglow model. 
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Fig. 5. Light Curve of the first 200 seconds of XRT observations. Labels mark the time 
segments in which the five spectra were extracted. 



analysis. We used for this analysis 5 XRT-WT spectra, selected as illustrated in figure 
El 2 XRT PC spectra, corresponding to the two Swift orbits where XRT operated in PC 
mode (see Fig.^, and PN and MOS spectra covering the full XMM- Newton observations. 
Results from XMM-Newton observations were also given by De Luca et al. (2005). 
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We first fitted the spectra with a simple power law combined with photoelectric 
absorption. We indicate with F = /3 + 1 the photon index. The Galactic column density 
along the line of sight to GRB 050713 A is TV^ 1-1 x 10^^ cm'^ (Dickey & Lockman 
1990). The results of our fits are shown in Table 2. Large spectral variations are evident, in 
particular between the spectra corresponding to the rise and fall of the first peak. Figure 
El shows the contours in the F — A^h plane for the five spectra in WT mode. The 
power law slope of the spectra WT2a and WT2b differ by AF « 1. It is also interesting 
to note that, during the decay WT2b, the spectrum recovers the same spectral slope as 
during the first decay WTl. This seems to indicate that the first decay WTl is not the 
soft X-ray counterpart of the prompt GRB emission but it is actually the tail of a fiare 
coincident with the flare detected by BAT in the 15-100 keV band 53 s after the Iq. 
In fact, as we see from Table 2, the spectrum of WTl is softer than the typical GRB 
spectrum. A gradual variation of the spectral power law index is also evident from WTl 
to WT4, where the spectral index is similar to that observed at later times by XRT in 
PC mode and by XMM-Newton. 

Motivated by the synchrotron emission model (Sari, Narayan & Piran 1998), we then 
fitted the spectra using a broken power law model. We kept the low energy spectral power 
law index fixed aX j3i — T i ~ 1 — 0.5, as expected from standard synchrotron models and 
left the high energy power law index /?2 = F2 — 1 as a free parameter. The results are 
shown in Table 3. The of spectra WTl, WT3, and WT4 are indistinguishable from 
those of Table 2. For spectra WT2a and WT2b the improvement in is significant at 
the 2% and 0.15% confidence level, respectively, using the F test. 

In both sets of fits there is a marginal evidence (between 2 and 3 a) that the absorbing 
column at 10** — 4 x 10"* seconds after the trigger was smaller than during the first 100-150 
seconds. 



Table 2. Single power law fits 



Spectrum 


Nh 
10^^ cm^^ 
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X'(dof) 


WTl 


0.61±0.10 


2.68±0.23 


51.9(45) 


WT2a 


0.53±0.09 


1.60±0.15 


76.4(69) 


WT2b 


0.62±0.05 


2.60±0.15 


156 (123) 


WT3 


0.44±0.08 


2.52±0.22 


53.1(42) 


WT4 


0.34±0.08 


2.09±0.24 


20.1(34) 


PCI 


0.30±0.16 


1.66±0.50 


5.9(4) 


PC2 


0.29±0.18 


2.09±0.40 


7.2(7) 


PC1+PC2 


0.28±0.08 


1.75±0.22 


31.9(14) 


XMM 


0.31±0.02 


2.06±0.04 


749(715) 
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Table 3. Broken power law fits 



Spectrum 
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The statistics of the PC spectra is not good enough to provide robust results adopting models 
more complex than a single power law. 
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Fig. 6. a) Left panel: single power law F — A^h confidence contours for the five WT 
spectra; b) right panel: broken power law r2 — inbreak confidence contours for the five 
WT spectra. 

5. Late-time "energy injection"? 

Rebrightenings and bumps in the afterglow light curve are generally due to either density 
inhomogeneities in the interstellar medium or to "energy injections" at later times. The 
latter could be produced by slower shells that catch up with the afterglow shock at later 
times. We will consider both possibilities and show that the energy injection model is 
favored by the data. Another possibility is that such flares are produced by late internal 
shocks, i.e. slower shells colliding with each other before deceleration (e.g. Fan & Wei 
2005; Zhang et al. 2005). 

First, let us consider the single powerlaw fit to the spectral data. During the time 
decay phases (time intervals WTl and WT2b in Figure EJ, the single power law spec- 
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trum (Table 2) is consistent with cx i/^^-^. On the other hand, during the rise phase 
WT2a, cx i/^"-^ as in the standard synchrotron model (e.g. Sari, Narayan & Piran 
1998). In this model, the changes of spectral slopes suggest that during the rise time 
the synchrotron frequency Vm sweeps across the observation band (0.5-10 keV). This is 
impossible to achieve with a density bump, since Vm is independent of the density for a 
completely adiabatic evolution, and almost independent of the density n [vm oc n^^^^^) in 
a fully radiative evolution. Therefore, we consider the possibility that the rebrightenings 
are the result of energy injections. 

A simple argument then shows that the synchrotron frequency is likely to vary within 
the 0.5-10 keV band covered by the XRT observation, and not be outside of the band 
as implied by the single powerlaw fit to the data. Indeed, for a completely adiabatic 
evolution, one has oc E^'^ (where E is the burst energy), while for a fully radiative 
evolution Vm oc £^*/^, both implying that the flux rebrightening F2/F1 scales with the 
frequency change factor as {I'm.^/vm,!)^ ■ If the synchrotron frequency were outside the 
observation band (i.e. at an energy ^0.5 keV before the beginning of the rebrightening, 
and Z 10 keV after the end of the rebrigthening) , then it would need to vary by a factor 
Ji2Q, implying an energy injection [vmal ^m,!)^ > 400. The rebrightening implied by 
this model is much larger than that seen in the data, F2/F1 ^ 10, implying that the 
synchrotron frequency lies within the 0.5-10 keV band. For this reason we consider more 
physical the broken power law fits in Table 3. The most striking feature from these fits 
is that the break frequency is smaller than 1.6 keV in the WTl spectrum and increases 
up to several keV in the WT2a spectrum (flux rebrightening) and then decreases again 
at energies < of 1-2 keV at later times. This is consistent with the flat optical to X-ray 
spectral index from Figure at 150-190 seconds after the trigger, indicating that, if the 
fast cooling regime applies, at this time the synchrotron frequency is between the X-ray 
and the optical band. 

For the same reasons discussed above, the shift of the synchrotron frequency could 
not have been caused by a density bump. On the other hand, in the energy injection 
scenario, a flux rebrightening by a factor of 10 would have been produced by an increase 
in break frequency by a factor of about 3, fully consistent with the data. 

Let us consider now the temporal behavior of the rebrightening as a function of the 
energy (sections WT2a and WT2b of the light curve in figureEjl. Our analysis shows that 
the temporal power law index during the rise time (section WT2a) is consistent with 
the constant value ~ 0.6 in all the spectral bands. During this phase, the temporal 
evolution is indeed dominated by the time-dependent energy supply, which is independent 
of the frequency band. The temporal flux decay index ad during the decay phase WT2b is 
shown in Figure|3| The lowest and highest energy bands are consistent with straddling the 
synchrotron frequency Vm for an index p — 2 of the power-law distribution of electrons, 
further confirming the broken power-law interpretation. Let us clarify this point. In the 
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low energy band (< 2 kev), the flux can be written as F^, oc (!^/i^c)^^^^^i/,max oc t"^ , 
while in the high energy band F^, cx {vm / Vc)^^^"^ {v / Vm)^^^'^ Fy^max oc i"^ where ai and 
a2 are the power law indices of the decay in the low and high energy band respectively 
(Sari, Piran & Narayan 1998). We get a good fit to the data with ai ~ —0.25 and 
a2 = — 3p/4 + 1/2. The dashed line in figure shows the results of this model for two 
values of p, p = 2 and p = 2.5. 

Taken at face value, the value of p which fits best the decay indices in figure 01 is 
somewhat different from the value implied by the F which fits best the spectrum during 
the decay phase WT2b, p ^ 3. However it should be considered that our description of 
the data is clearly over-simplified, and that statistical and systematic uncertainties in 
our determination of and F can certainly alleviate this discrepancy. In conclusion, we 
find remarkable that, at least qualitatively, the energy injection model in the framework 
of the afterglow theory is able to reproduce both the spectral and the temporal behavior 
observed in GRB 050713A. 

The time lags between the light curves in the six energy bands considered here can be 
further used to constrain the product of the magnetic field energy fraction e^, the burst 
energy, E, and the density of the external medium, n, e^'^i^gj^n^^, using the expression 
of the synchrotron cooling time given by Sari, Piran & Narayan (1998). The two curves 
in Figure ^ show the expected time lags for e^'^-Eg^^n^'^ = lO^'^ and 2 x 10^^. It is 
interesting to note that the synchrotron model does predict a time lag of the order of 
6-9 s between the low energy photons (< 1 keV) and higher energy photons (2-10 keV), 
and a flattening of the time lag for energies > 2 keV as observed. 

6. Summary 

We have presented time and spectral analysis of GRB 050713A, whose light curve dis- 
plays several rebrigthcnings both in the 7-ray and in the X-ray bands. Our time-resolved 
spectral analysis has allowed us to conclude that the rebrigthcnings are most likely due to 
late energy "injections" . Since the main energy output of the GRB engine rapidly decays 
with time (Janiuk et al. 2004), the energy injections are likely due to either slower shells 
generated during the prompt phase of the GRB engine and catching up at later times, 
or to later energy production by, for example, disk fragments that accrete at later times 
(Perna et al. 2006). 

We would like to emphasize two conclusions that came out from the analysis of the 
light curve of GRB050713A. First the behavior of the time lags between the light curves in 
the six energy bands is consistent with what expected from synchrotron cooling. Second, 
the spectral index of the first decay spectrum WTl is consistent with that of the decay 
from the main fiare WT2b, thus suggesting that the WTl is actually the tail of a flare 
coincident with the flare detected by BAT. 
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The X-ray light curve may present features detected in many Swift bursts (Chincarini 
et al. 2005, Nousek et al. 2005, Burrows et al. 2005b). In particular we notice that the 
XRT decay index starting 200 s after the trigger is flatter than the decay index during the 
XMM observation at about 8 hours from the trigger. This implies that a break should 
occur between 5000 s and 20.000 s after the trigger. The decay index after the break 
appears consistent with the predictions of the uniform ISM afterglow model, while the 
decay index before the break is much shallower. This shallow-to-normal decay may be 
due to the cessation of the refreshed shock phase. However we cannot exclude that this 
is a jet break due to the deceleration of the outflow. The break occurs at a time when 
the bulk Lorentz factor of the shock has slowed to F ^ 1/6*^6* (Rhoads 1997). 

All optical points seem to be typical of the afterglow emission produced by an external 
forward shock. The optical flux varies with time as F^, ^ which is consistent with 

the afterglow decay during the radiative evolution, Fj, ~ t^^^'^ . We should notice that 
this temporal behavior is similar to the X-ray one between WT4 and XRT-PC. 
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